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Abstract 

Electrolytic capacitors are used in several applications rang- 
ing from power supplies for safety critical avionics equipment 
to power drivers for electro-mechanical actuators. Past expe- 
riences show that capacitors tend to degrade and fail faster 
when subjected to high electrical or thermal stress condi- 
tions during operations. This makes them good candidates for 
prognostics and health management. Model-based prognos- 
tics captures system knowledge in the form of physics-based 
models of components in order to obtain accurate predictions 
of end of life based on their current state of health and their 
anticipated future use and operational conditions. The focus 
of this paper is on deriving first principles degradation mod- 
els for thermal stress conditions and implementing Bayesian 
framework for making remaining useful life predictions. Data 
collected from simultaneous experiments are used to validate 
the models. Our overall goal is to derive accurate models of 
capacitor degradation, and use them to remaining useful life 
in DC-DC converters. 

1. Introduction 

Avionics systems play a critical role in many aspects of air- 
craft flight control. As the system complexity and flight crit- 
icality of functions performed by these systems increases, 
the related consequences of in-flight malfunctions are bound 
to increase. This drives the need for Integrated Vehicle 
Health Management (IVHM) technologies for flight-critical 
avionics. Studying and analyzing the performance degrada- 
tion of embedded electronics in the aircraft domain is abso- 
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lutely necessary to increase aircraft reliability, assure in-flight 
performance, and reduce maintenance costs, (J. R. Celaya, 
Wysocki, Vashchenko, Saha, & Goebel, 2010; Ferrell, 1999). 
In addition to this, an understanding of the behavior of de- 
teriorated components is needed as well as the capability to 
anticipate failures and predict the remaining useful life (RUL) 
of the electronic systems. 

An avionics system module consists of hardware (power sup- 
ply, Global positioning system (GPS) receiver, Intertial mea- 
surement unit (IMU)) and software (GPS software, INAV - in- 
tegrated navigation solution) components (Kulkami, Biswas, 
Bharadwaj, & Kim, 2010). Switched-mode power supplies 
are widely used in DC-DC converters because of their high 
efficiency and compact size. Buck-boost DC-DC converter 
steps voltage levels up/down based on the application re- 
quirements, in which electrolytic capacitors and metal ox- 
ide semiconductor field effect transistor (MOSFET) switches 
are known to have the highest degradation and failure rates 
among all of the components (Goodman, 2005). Degraded 
capacitors affect the performance and efficiency of the DC- 
DC converters in a significant manner and also impose a risk 
on instantiating cascading failures on other connected subsys- 
tems. 

In this paper we develop an effective prognostics and health 
management (PHM) methodology that applies to electronic 
systems and components. In particular, we develop a method- 
ology to enable early detection of failure precursors in capaci- 
tor elements associated with DC-DC power supplies. Our ap- 
proach combines physics-based modeling supported by em- 
pirical experimental analysis for improving the degradation 
models, and then implementing bayseain framework to pre- 
dict remaining useful life of electrolytic capacitors. Our hy- 
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pothesis is that early detection will lead to better prediction 
and end of life estimates by tracking and modeling the degra- 
dation process. 

The structure of the paper is as follows. Section II dis- 
cusses the prognostics metholdology implemented for this 
work. Section III presents introduction to electrolytic capaci- 
tors and its basic structure, operation and degradation mech- 
anisms. Section IV discusses capacitor first principle models 
in detail. Section V describes the thermal stress aging exper- 
iments conducted for this work. Section VI and VII presents 
the prognostic framework methodology and RUL results re- 
spectively. The paper ends with discussion and conclusion in 
section VIII. 

2. Prognostic Methodology 

Prognostics is the process of predicting health condition and 
remaining useful life based on current and previous state, cur- 
rent and future operating conditions. Prognostics and health 
management (PHM) methods combine sensing, data collec- 
tion, interpretation of environmental, operational, and perfor- 
mance related parameters to indicate systems health as well 
as anticipate damage propagation due to degradation. PHM 
methodologies can be implemented through the use of var- 
ious techniques that study parameter variations, which indi- 
cate changes in performance degradation based on usage du- 
ration and conditions. 

Prognostics is an essential technology for improving system 
safety, reliability, and availability. Prognostics deals with de- 
termining the health state of components, and projecting its 
health evolution into the future to make end of life (EOL) 
and remaining useful life (RUL) estimations. Model-based 
prognostics approaches perform these tasks with the help of 
first principles based physics models that capture knowledge 
about the system, its components, and their failures (Daigle 
& Goebel, 2011; Saha & Goebel, 2009). 

We adopt an approach wherein we develop detailed physics- 
based models of components and systems that include de- 
scriptions of how fault parameters evolve in time. The im- 
plemented prognostics architecture is as shown in Figure 1 . 
We have implemented this approach on empirical degradation 
models in our earlier work (J. Celaya et al., 2011a; Kulka- 
mi et al., 2012). Experimental studies are being conducted 
and first principles based degradation models are being de- 
rived using the descriptions mentioned in (Kulkami, Celaya, 
et al., 2011; J. Celaya et al., 201 la; Fife, 2006; MIL-C-62F, 
2008). Identifying the failure precursors and developing ac- 
curate models of degradation/ failure has been the most diffi- 
cult problem of our research goal. These models depend on 
known as well as unknown and possibly time-varying wear 
parameters. Early detection and analysis may lead to better 
prediction and end of life estimates of the capacitor by track- 
ing and modeling the degradation process. Faults and degra- 


dations appear as parameter value changes in the model, and 
this provides the mechanisms for tracking system behavior 
under degraded conditions (J. Celaya et al., 201 la, 201 lb). 
The derived state space models are then implemented in a 
Bayesian framework for prognostics. 



Figure 1 . Prognostics Methodology 


In the next section we discuss in brief the basics of elec- 
trolytic capacitors, their detailed structure and different 
mechanisms under which the devices degrade. 

3. Electrolytic Capacitors 

An aluminum electrolytic capacitor, as illustrated in Figure 2, 
consists of a cathode aluminum foil, electrolytic paper, elec- 
trolyte, and an aluminum oxide layer on the anode foil sur- 
face, which acts as the dielectric. When in contact with the 
electrolyte, the oxide layer possesses an excellent forward di- 
rection insulation property (Bengt, 1995). Figure 3 shows a 
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detailed view of the cross section of an electrolytic capacitor 
structure. To get higher capacitance values for the same sur- 
face area of the anode and cathode foils, the foil is etched by a 
chemical process. Together with magnified effective surface 
area attained by etching the foil, a high capacitance value is 
obtained in a small volume (Fife, 2006). Since the oxide 
layer has rectifying properties, a capacitor has polarity. If 
both the anode and cathode foils have an oxide layer, the ca- 
pacitors would be bipolar. In this work, we analyze non-solid 
aluminum electrolytic capacitors in which the electrolytic pa- 
per is impregnated with liquid electrolyte. After etching, the 
plates are anodized by coating them with a thin aluminum ox- 
ide layer on the surface of the foil. This layer of aluminum 
oxide acts as the dielectric (insulator) and serves to block the 
flow of direct current between the two capacitor plates (Fife, 
2006). 


of functionality due to a short or open circuit, and (2) degrada- 
tion failures, where there is gradual deterioration of capacitor 
due to accumulated damages. 

The fishbone diagram in Figure 6 summarizes the most com- 
mon set of failure modes for electrolytic capacitors that have 
been discussed in (J. Celaya et al., 2011a; Kulkami et al., 
2012). This diagram identifies the relationship between root 
causes and failure modes observed in electrolytic capacitors. 
These root causes can occur individually or combined man- 
ner in a capacitor depending upon the conditions of operation. 
Our focus in this work is on the thermal stressors that govern 
the capacitor degradation, specifically, we study high temper- 
ature scenarios and their effects on the electrolytic capacitor 
degradation. 

3.1. Equivalent Electrical Circuits 
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Figure 3. Capacitor Detail Structure 


Electrolytic capacitor performance is strongly affected by its 
operating conditions, such as voltage, current, frequency, and 
ambient temperatures. Degradation in the capacitor manifests 
an increase in the equivalent series resistance ( ESR ) and de- 
crease in capacitance (C), due to deterioration of electrolyte 
quality, decreases in electrolyte volume due to evaporation, 
weakening of the oxide layer, over operating time (Bengt, 
1995; Fife, 2006). The ESR of a capacitor is the sum of 
the resistance due to aluminum oxide, electrolyte, spacer, and 
electrodes (foil, tabbing, leads, and ohmic contacts) (Bengt, 
1995) and capacitance is the ability of a capacitor to store 
charge in an electric field. The health of a capacitor is often 
measured by the values of these two parameters. 

There are certain industry standard thresholds for these pa- 
rameter values, if the measurements exceed these thresholds 
then the component is considered unhealthy, i.e., the compo- 
nent has reached its end of life, and should be immediately 
replaced before further operations (Lahyani, Venet, Grellet, 
& Viverge, 1998). Failures in a capacitor can be one of two 
types: (1) catastrophic failures, where there is complete loss 


A simplified electrical lumped parameter model of 
impedance, M\, defined for an electrolytic capacitor is 
shown in Figure 4. The ESR dissipates some of the stored 
energy in the capacitor. An ideal capacitor would offer no 
resistance to the flow of current at its leads. 

11 ESR 
•— 1| — WW 

C 

Figure 4. Lumped Parameter Model { M\) 

It has been observed that under thermal overstress storage 
conditions (Bengt, 1995; J. Celaya et al., 2011a), the capac- 
itance (C) and ESR value depends of the electrolyte resis- 
tance Re- A more detailed lumped parameter model derived 
for an electrolytic capacitor under thermal overstress condi- 
tion, M 2 can be modified from M\, as shown in Figure 5. Ri 
is the combined series and parallel resistances in the model. 
Re is the electrolyte resistance. The combined resistance of 
Ri and Re is the ESR of the capacitor. C is the total capac- 
itance of the capacitor as discussed earlier. 

ESR 

1 — ^ — 1 j. 

►-Wr AAA — — • 

Ri Re q 

Figure 5. Lumped Parameter Model {M 2 ) 

4. Physics Based Modeling 

In this section we discuss about deriving the first principles 
based degradation models for capacitors under thermal over- 
stress conditions. Under thermal overstress conditions since 
the device was subjected to only high temperature with no 
charge applied we observe degradation only due to electrolyte 
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Causes\ Failure Mechanisms 


Failure Mode 


Figure 6. Fishbone diagram of failure mechanisms in aluminum electrolytic capacitor 


evaporation. The degradation models are derived based on 
the underlying physics of operation and measurements from 
experimental data. 

4.1. Structural Model 

For deriving the physics based models of an electrolytic 
capacitor it is also necessary to know about the structural 
and manufacturing details, since health estimations are done 
based on the type of electrolyte, volume of electrolyte, oxide 
layer thickness etc. The models defined use this information 
for making effective degradation/failure predictions. A de- 
tail structural study of the electrolytic capacitor under test is 
discussed in this section. 

During modeling it is not possible to know the exact amount 
of electrolyte present in a capacitor. But using information 
from the structural details as shown in Figure 7, we can ap- 
proximately calculate the amount of electrolyte present. A 
very highly porous separator paper is used which soaks all 
the electrolyte. The paper separator is sandwiched between 
anode and cathode, each having a thickness ds , d A and dc 
respectively ( ds ~ dA ~ dc ). Based on the type and configu- 
ration, the electrolyte volume will vary which can be updated 
in the model parameters. 

The equation for calculating the approximate electrolyte vol- 
ume is derived from calculating the total capacitor capsule 
volume, V c given by : 

V c = ir r\h c ( 1 ) 

where : 

r c = radius of capacitor capsule 
h c = height of capacitor capsule 



The approximate electrolyte volume, V e based on all the other 
known structural details of the capacitor can expressed as: 

Ve = 7T r 2 c h c - A s (d A 4- d c ) (2) 

4.2. Capacitance Degradation Model 

Exposure to high temperatures, T avp iied > T ra ted results in 
accelerated aging of capacitors (Kulkami, Celaya, et al., 
2011; J. Celaya et al., 2011a; 60068-1, 1988). Higher 
ambient storage temperature accelerates the rate of elec- 
trolyte evaporation leading to degradation of the capaci- 
tance (Kulkami, Celaya, et al., 2011; Bengt, 1995). The de- 
pletion in electrolyte volume, V e , (Kulkami, Biswas, et al., 
2011; Rusdi et al., 2005) is given by : 

v e (t) = V 0 - (A s j e0 w e )xt (3) 

where: 

V 0 = initial electrolyte volume 
j eo = evaporation rate (mg ram -1 area~ l ) 
w e = volume of ethyl glycol molecule 
t = time in hours. 

The total lumped capacitance from first principles of electro- 
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magnetism is given by : 


From Eq. (9) we can express Eq. ( 8 ) as : 


c = ( 2e R e 0 A s )/dc (4) 


where: 

6 r = relative dielectric constant 
cq = permitivity of free space 


^ ^ dC A 

CWi a — 4- -j— A t, 

at 

C k+ \a = C k a - ( w e A s j eo )At , hence 

C k+1 =C k - {We ^ Jeo) At. 

a 


( 10 ) 


Thus from Eq. (3) and Eq. (4) we can derive the first princi- 
ples capacitance degradation model, M 3 given by : 


M 3 : C(t ) = 


f 2efieo\ 
V dc ) 


Vo-Ve{t) \ 

Jeo t W e ) 


(5) 


As discussed earlier, increase in the core temperature evapo- 
rates the electrolyte thus decreasing the electrolyte volume 
leading to degradation in capacitance. The resultant de- 
crease in the capacitance can be computed using model, M 3 
wherein the decrease in electrolyte volume, (Ve) leads to de- 
crease in capacitance, (C). 


4.3. Dynamic Models 


The complete discrete time dynamic model for capacitance 
degradation can be summarized as : 

M 4 : C k+1 =C k - W eA sJeo) At 

a c 

Model A^ 4 , in Eq. (1 1) is implemented in a Bayesian track- 
ing framework. In this work we are implementing a un- 
scented Kalman filter (UKF) since the degradation in capac- 
itance (state) due to decrease in electrolyte is considered to 
be a dynamic linear model and the evaporation rate (J eo ) pa- 
rameter, assumed to be varying and estimated online. Next 
we discuss the implementation of the Bayesian framework 
methodology for prognostics (Saha & Goebel, 2009; Daigle 
et al., 2012; Daigle & Goebel, 201 1 ). 


The electrolyte volume. 14 can also be calculated and can also 
be expressed as V e = A s .dc . Hence the oxide surface area, 
A s can be expressed in terms of electrolyte volume, V e and 
oxide thickness, dc as: 

1/ 

A s = (6) 


From Eq. (4) and Eq. ( 6 ), dynamic capacitor degradation 
model can be updated as : 




( 2 e R e 0 

V dc 



(7) 


From Eq. (3) the first order discrete approximation for change 
in electrolyte volume can be expressed as: 


^ — ( w e A s j eo ), 

dVe 

Ve(k+ 1 ) = V e (k) + £, 

Ve(/c+l) Ve(k) (w e A s j eo )At. 


( 8 ) 


From Eq. (7) we have, 


5. Unscented Kalman Filter 

Estimation in nonlinear system is very important because 
many practical systems involve nonlinearities in their oper- 
ation through one form or another. Estimation of the state 
accurately of such non-linear system is very important to di- 
agnostics and be further implemented to prognostic appli- 
cations which is extremely difficult. The Extended Kalman 
Filter (EKF) which applies the KF to nonlinear system by 
linearizing all nonlinear models, has become a most widely 
used method for estimation of nonlinear system. Although 
the EKF maintains is computationally efficient recursive up- 
date form of the KF for non-linear systems, it suffers a num- 
ber of serious limitations (Julier & Uhlmann, 1997, 2004). 

1 . Only reliable if the error propagation to the future states 
can be approximated by a linear function. 

2. Linearization can be applied only if the Jacobian matrix 
exists i.e, cannot be implemented if the system has dis- 
continuities. 

3. Calculating Jacobian matrices can be a very difficult and 
error-prone process. 


where: 


V e(k) 

V e(k) 


(C k )a. 


2 

c> 


d 2 c 

2e/?eo 


The Unscented Kalman Filter (UKF) was proposed by (Julier 
& Uhlmann, 1997, 2004) to overcome these problems in non- 
linear systems. The unscented Kalman filter, instead of ap- 
proximating the nonlinearity, approximates the state distribu- 
tion (Julier & Uhlmann, 1997, 2004). This procedure main- 
tains the nonlinear functions exactly, eliminating the need to 
calculate Jacobian’s, and thereby providing an easier imple- 
mentation framework. In this section we will look in the 
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detail framework of UKF which has been implemented for 
prognostics in this work. 

A nonlinear system, described by the difference equation and 
the observation model with additive noise is given as : 

Xk=~P[xk-i,Uk-i,k\+Wk-i ^ 

Zk = H[xA;,fc] +Vk 

where x(k) is the n-dimensional state of the system at time 
step k , u(k) is the input vector, w(k) is the process noise, 
Q, z (k) is the observation vector and Vk is the measurement 
noise, R. 

To solve the problem of predicting the future state or obser- 
vation based on the Unscented transform (UT), UT takes a 
random variable x , with mean x and covariance J? xx , which 

is related to a second random variable y by a nonlinear func- 
tion y = f{x). A small set of points ( sigma points) with 
mean y and covariance P TX are selected (Julier & Uhlmann, 
1997), which are deterministically selected and weighted to 
exhibit properties to match the mean and covarinace of the 
original distribution. A non-linear transformation is applied 
to each point to get the transformed points, statistics of the 
transformed points is then calculated to estimate the mean 
and covariance of the transformed points. The sigma point 
weights do not directly represent probabilities of the sigma 
points, and hence do not have to lie in the interval [0, 1). The 
weights Wi can be positive or negative, but need to obey the 
following conditon to provide an unbiased estiamte. 

Yh Wx = l (13) 


Each sigma point is instantiated through the function(/) to 
obtain new set of sigma points 3^- 

= f{Xi) (i4) 


Uhlmann, 1997, 2004). At each step, the unscented trans- 
form is applied to the state estimate and is used for a single 
step prediction. In contrast, here, we apply the transform to 
the state parameter distribution at given single time point kp 
, and use this for multi-step predictions to EOL. There are 
several methods which exits for selecting the sigma points 
out of which we implement the symmetric unscented trans- 
form for the prognosis problem (Daigle et al., 2012; Julier & 
Uhlmann, 2004). Detailed results will be presented in Section 
6 for capacitor degradation problem. 

In the symmetric unscented transform, 2 n x + 1 sigma points 
are selected symmetrically about the mean as follows (Julier 
& Uhlmann, 2004): 


Ab 

Wo 

Wi 


*i + n 

W i+n 


n + k 

x + \J ((n -f- k)P xx )i, 
k 

2 (n + k ) 

X y/((n “I - k)P xx )i 
k 

2(ti ~b k) 


(17) 


where ^/((n + k)P xx )i refers to the i th column of the ma- 
trix square root of (n + k)P XXl computed using the Cholesky 
decomposition since it is numerically efficient and stable. Pa- 
rameter k is used to tune the higher moments of distribution 
and suggested to have a smaller values as possible to bring the 
sigma points closer together. If x is assumed Gaussian, then 
selecting n + k = 3 is recommended (Julier and Uhlmann, 
1997). 

6. Thermal Overstress Experiment 


The mean of the transformed points is given by: 

v = 05 ) 

i 

The covariance of the transformed points is given by: 

Pyy “ W>)(*0 -M°) T (16) 

i 

The basic idea of the unscented transform is that it is eas- 
ier to apporximate a probability distribution x than it is to 
approximate an arbitary nonlinear function / or transforma- 
tion (Julier & Uhlmann, 2004). This basic principle is imple- 
mented in the unscented Kalman filter where the unscented 
transform is exploited for nonlinear state estimation (Julier & 


In this setup we emulated conditions similar to high tempera- 
ture storage conditions (Kulkami, Celaya, et al., 2011; Kulka- 
mi, Biswas, et al., 2011), where capacitors were placed in 
a controlled chamber and the temperature raised above their 
rated specification (60068-1, 1988). Pristine capacitors were 
taken from the same lot rated for 10V and maximum storage 
temperature rating of 85°C. 

Experiments were conducted with 2200 pF capacitors with 
TOS temperature at 105°C and humidity factor at 3.4%. 
The chamber temperature was gradually increased in steps of 
25°C till the pre-determined temperature limit was reached. 
The capacitors were allowed to settle at a set temperature for 
15 min and then the next step increase was applied. This 
process was continued till the required temperature limit was 
attained. To decrease possibility of shocks due to sudden de- 
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crease in the temperature the above procedure was followed. 
At the end of specific time interval the temperature was low- 
ered in steps of 25° C till the required room temperature was 
reached. 

Before being characterized the capacitors were kept at room 
temperature for 15 min. The ESR value is the real impedance 
measured through the terminal software of the instrument. 
Similarly the capacitance value is computed from the imag- 
inary impedance using Electrochemical Impedance Spec- 
troscopy (EIS). Characterization of all the capacitors was 
done for measuring the impedance values using an SP- 
150 Biologic impedance measurement instrument (Biologic, 
2010). Figure 8 shows the plots decrease in capacitance due 
to accelerated aging for all the 15 capacitors under test at dif- 
ferent aging times. 
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(b) Output error for Cap # 5 
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Figure 8. Capacitance Plot for all the devices under TOS 


In the thermal overstress experiments, the capacitors we char- 
acterized periodically and after 3400 hours of operation it was 
observed that the average capacitance ( C ) value decreased by 
more than 9-11% while decrease in ESR value was observed 
around 20 - 22%. From literature (60068-1, 1988) under ther- 
mal overstress conditons higher capaitance degradation is ob- 
served and minor degradation in ESR which correlated with 
the data collected. The failure thresholds under storage condi- 
tions for capacitance (C) is 10% while that for ESR is around 
280- 300% of the pristene condition values (60384-4-1, 2007; 
Kulkarni, Celaya, et al., 2 01 1). Hence the capacitance degra- 
dation data was used as a precursor to failure parameter to 
estimate the current health condition of the device. 

7. Prediction of Remaining Useful Life Results 

State estimation and RUL estimation results are discussed for 
capacitor Cap # 5 out of a batch of 15 available capacitors 
under test. Figure 9 shows the result of filter tracking for 
degradation in capacitance upto 3200 hours of aging time. 


Figure 1 0 presents results from the remaining useful life pre- 
diction algorithm at different aging times t p = 87, 607, 
1495, 2131,2800 (hrs), at which the capacitors are charac- 
terized and their capacitance (C) value is calculated. The 
failure threshold is considered to be 10% decrease in ca- 
pacitance value, which in this case is at 3200 hours of ag- 
ing time. End of life (EOL) is defined as the time at which 
the forecasted capacitance value trajectory crosses the EOL 
threshold. Therefore, RUL is EOL minus aging times t v = 
87.5, 607, 1495, 2131, 2800 (hrs). 
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An a-X prognostics performance metric (Saxena et al., 2009, 
2008) is presented in Figure 1 1 for test case of Cap #5. The 
central dashed line represents ground truth and the shaded 
region is corresponding to a 30% (a = 0.3) error bound in the 
RUL prediction. Performance metric identifies whether the 
algorithm performs within desired error margins (specified by 
the parameter a) of the actual RUL at any given time instant 
(specified by the parameter A) (Saxena et al., 2009) and is 
based on relative accuracy (RA) metric in Eq. (18). 


RA = 100 


RUL* -RUL'\ 

RUL* J 


(IB) 


Table 1 shows the performance summary based on the RA of 
all the capacitors under thermal stress performance. These 
metrics allows for an assessment of the percentage accuracy 
relative to the ground-truth value. RA values of 100 represent 
perfect accuracy. The RA is presented for all the test cases 
for different prediction times. The last column of Table 1 
represents the median RA of all the test cases for a particular 
prediction time. It must be noted that if the prediction error 
magnitude grows beyond 100% RA gives a negative value. 
We do not consider such cases since these cases would not 
have qualify the tests for calculating RA (Saxena et al., 2009), 
these are indicated by NA in Table 1 . 

From the a-X metric plot in Figure 1 1 it can be observed that 
the relative accuracy is not as good at the end but the accuracy 
is good enough under acceptable limits. This is due to the 
non-linearity observed in the data at the end of the aging time 
and the limitation of the model due to not including the oxide 
layer breakdown. The residuals show an increased error with 
aging time, since the breakdown in the oxide layer observed 
due to stress is not considered for this model which starts 
to dominate in the later stages of aging of the device. This 
breakdown is exponential in nature and as we can observe a 
dip in the capacitance values from the linear path in the later 
stages. 

8. Conclusion and Discussion 

This paper presents a first principles based degradation elec- 
trolytic capacitor model and an parameter estimation algo- 
rithm to validate the derived model, based on the experimen- 
tal data. The majors contributions of the work presented in 
this paper are: 

1. Development of the first principles degradation model 
based on accelerated life test aging data which includes 
decrease in capacitance as a function of time and evapo- 
ration rate linked to temperature conditions; 

2. Implementation of a Bayesian based health state track- 
ing and remaining useful life prediction (RUL) algorithm 
based on the UKF filtering framework; 

3. Prediction of remaining useful life for capacitors based 
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Figure 11. Performance based on Alpha-Lambda metric for 
Cap#5 


first principles degradation model, .M3; 

The degradation model, M 3 based on the first principles 
gives an indication of how a specific device degrades based 
on its structure, material properties, operating conditions, etc. 
The results presented here are based on accelerated aging ex- 
perimental data and on the accelerated life timescale. In our 
earlier work we studied the degradation models based on the 
observed data, and the work discussed here is a next step to 
generalize the model. Though as discussed in section 4, as a 
first step a dynamic linear model has been implemented for 
degradation model. This degradation model for decrease in 
capacitance, C and varying evaporation rate j eo needs to be 
updated and include the model of break-down in the oxide 
layer which is exponential in nature and dominates in the 
later stages of aging. Further research will focus on devel- 
opment of functional mappings that will translate the acceler- 
ated life timescale into real usage conditions timescale, where 
the degradation process dynamics will be slower, and sub- 
jected to varying stress conditions. 

The performance of the proposed first principles degradation 
model, M 3 is acceptable for the current study based on the 
quality of the model fit to the experimental data and the RUL 
prediction performance from a-X metric plot. Additional ex- 
periments are currently underway to increase the number of 
test samples. This will greatly enhance the quality of the 
model, and guide the exploration of additional degradation- 
models, where the loading conditions and the environmen- 
tal conditions are also accounted for towards degradation dy- 
namics. 
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Nomenclature 

6r relative dielectric constant 

to permitivity of free space 

V Q initial electrolyte volume 

j eo evaporation rate (mg min~ l area ~ l ) 
p E electrolyte resistivity 

A s effective oxide surface area 

w e volume of ethyl glycol molecule 

d,A thickness of anode strip, 

dc thickness of cathode strip 

d s thickness of paper spacer 

M i electrical lumped parameter model 

M 2 updated lumped parameter model 
M$ capacitance degradation model 

M\ capacitance discrete time model 
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Tab le 1. Summary of RUL forecasting results 

RA a is the mean relative accuracy of all capacitors at each prediction time ( t p ) 
RAb is the mean relative accuracy of each capacitor at all prediciton times 
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